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Elevated temperature deformation processing, “hot working," is central to the pro- 
duction of more than 80% of all metal products. The advent of large deformation finite 
element methodologies is beginning to permit the numerical simulation of hot working pro- 
cesses whose design until recently has been based on prior industrial experience. Proper 
application of such finite element techniques requires realistic constitutive equations which 
more accurately model material behavior during hot working. 

The constitutive equations should satisfy several requirements. First, the equations 
should be able to model large, three dimensional deformations. Second, the model should 
encompass the range of hot working conditions, which includes strain rates ranging from 
10“ s to 10 1 sec -1 or greater, homologous temperatures from .5 to .9, and interrupted de- 
formation histories. Third, determination of parameters associated with the model should 
be straightforward and require a minimum of materials testing. Fourth, the model should 
provide a means of representing material microstructural state and its evolution during 
deformation. Finally, the model should be formulated with due consideration of issues 
regarding their numerical implementation in finite element programs. 

A simple constitutive model for hot working which satisfies most of these require- 
ments is the single-scalar internal variable model for isotropic thermo-elasto-viscoplasticity 
proposed by Anand [1985,1982]. In the next section we recall this constitutive model. The 
specific scalar functions for the equivalent plastic strain rate and the evolution equation for 
the internal variable presented here are slight modifications of those proposed by Anand 
[1982]. These modified functions are better able to represent high temperature material 
behavior. Following this presentation of the constitutive equations we briefly describe 
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our monotonic constant true strain rate and strain rate jump compression experiments 
on a 2% silicon iron. The material parameters appearing in the constitutive model can 
be determined from the stress-strain curves resulting from such experiments. The model 
is implemented in the general purpose finite element program ABAQUS [Hibbitt, et.al., 
1984]. Finally, using this program, the predictive capabilities of the model are evaluated 
for some simple deformation histories. 

Constitutive Model 


(a) Stress-strain-temperature rate relation: 
T* = .C[D-DP]-IW, 


( 1 ) 


where with T denoting the Cauchy stress and F denoting the deformation gradient, 


T = (det F)T 
fv eT-WT + TW 
£ = 2 mJ + [* ~ (2/3)m]1 ® 1 
m = A (*)» K = k(e) 

II = (3ica)l 
a 

D = sym FF _1 
W = skewFF -1 
0 


Kirchhoff stress; 

Jaumann derivative of Kirchhoff stress; 
Elasticity tensor; 

Shear and bulk moduli; 

Stress temperature tensor; 

Coefficient of thermal expansion; 
Stretching tensor; 

Spin tensor; 

Absolute temperature. 


(b) Flow rule: 


The constitutive equation for D p is: 

D’ = *(§£)■ 


where 

f — >0, o < s 

b = y/(3/2) T'-T' 


equivalent plastic tensile strain rate, 
equivalent tensile stress, 


( 2 ) 
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and s is a scalar internal variable with the dimension of stress. We call it the de- 
formation resistance which may be associated with a microstmctural characteristic 
such as dislocation density. 

(c) Evolution equation: 

The deformation resistance s is assumed to evolve according to: 
s = h(b t 0 t s)'f 

where h represents strain hardening (with dynamic recovery) and r represents static 
recovery. 

The major task in completing the constitutive model involves a specification of the 
rate equation / above and the evolution equation for the internal variable s. We propose 
the following model: 


f = Aexp (~S) [sinh({f)] 1/m , b < s, 

(3) 

« = ho[< 1 - (s/s*) >] a , a > 1, 

(4) 

with 


s* = s[£exp (ft)]\ 

(5) 


In this model the material parameters are A,Q , (,m,ho,a,s,and n, and k is Boltzmann’s 
constant. 

The static recovery function f in this model is set equal to zero. Accordingly, these 
equations are unable to model recovery during hold periods between high rate deformation 
passes or during very slow deformation processing. 

Compression Experiments on a 2% Silicon Iron 

Isothermal hot-compression tests have been performed to evaluate the material 
constants associated in the above model. A 2% silicon iron was selected since it does not 
dynamically recrystallize and since it does not experience a change in crystal structure 
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upon quenching from elevated temperature deformation. Tests were performed on a high 
temperature materials test system consisting of a vacuum furnace mounted in a servo- 
hydraulic test machine. An analog function generator was used to produce constant true 
strain rates and strain rate jumps. Data was collected on a microcomputer data acquisition 
system. 

The test specimens were circular cylinders with a height-to-diameter ratio 1.5. Shal- 
low, concentric, circular grooves were machined into the ends of each specimen to hold a 
high temperature lubricant. Lubricants were mixtures of powdered glass and boron nitride 
powder. Homogeneous compression to a true strain of approximately -1.0 was achievable 
on a routine basis. 

Two sets of experiments were performed: isothermal, constant true strain rate tests 
and strain rate jump tests. The constant true strain rate tests were performed for a 
range of temperatures from 800 to 1200 degrees Celsius, and a range of strain rates from 
10“* to 10° per second. Figure 1 provides a representative subset of these tests for a given 
temperature at varying true strain rates. Strain rate jump tests were performed to provide 
means of evaluating strain rate dependence at a given internal state , because such a test 
instantaneously decouples the strain rate equation (3) from the evolution equation for the 
interna] variable (4). Figure 2 shows a representative series of jump tests at a constant 
temperature. The stress/strain data following the jump in strain rate also provides data 
for an independent comparison of the constitutive model predictions and actual material 
response. 

Material Constant Determination and Model Evaluation 

Material constants for the model represented by equations (3) to (5) have been 
determined using data obtained from the constant true strain rate and strain rate jump 
tests. The equations have been incorporated using procedures outlined by An and [1986] 
via a user-material interface in the finite element program ABAQUS. Stress-strain curves 
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have been calculated using ABAQUS for conditions representative of the experiments on 
the silicon iron. The calculated curves are compared with the experimental data in Figures 
3 and 4. The agreement between the theory and experiment is excellent. 

Several test specimen geometries have also been evaluated which produce a gradient 
of internal microstructure within a single deformed specimen. These specimens have been 
deformed at temperature and then quenched. Different measures of microstructure, such as 
microhardness and etchpit density, have been compared with the variation in the internal 
variable predicted by the constitutive model. 
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FIGURE 4. Predicted Strain Rate Jump Test, 
Stress/Strain Curves versus 
Experimental Data . 
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